Invariant color calculus and generalized Balitsky-Kovchegov hierarchy 
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We derive generalization of the Balitsky-Kovchegov (BK) equation for a dipole, which consists 
of a parton and an antiparton of arbitrary charge. At first, we develop one method of indexless 
transformation of color expressions. The method is based on an evaluation of the Casimir operator 
on a tensor product. From the JIMWLK equation we derive the evolution equation for a single 
parton and prove gluon Reggeization in an arbitrary color channel. We show that there is a color 
duplication of such Regge poles. Higher t-channel color exchange has its own Regge pole, which 
residue is proportional to the quadratic Casimir. Taking a fundamental representation, we derive 
the usual BK equation and shed new light on the meaning of linear and nonlinear terms. Finally, 
we discuss a linearized version of the generalized BK equation. 

PACS numbers: 12.38.-t, 25.75.-q 



I. INTRODUCTION 



The Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner(JIMWLK) equation [3, Q is an important part of 
our understanding of high energy evolution of QCD scattering amplitudes. It describes a scattering of a dilute 
• projectile on a dense target. The main disadvantage is its functional form, which is necessary for studying of an 
Q-f arbitrary projectile. The functional equation is difficult to solve, even numerically. The Balitsky-Kovchegov (BK) Q 

■ hierarchy is a special case of the JIMWLK equation where the initial projectile is fixed and taken by a quark- antiquark 
pair. To perform practical calculations one can use the mean field approximation, which allows one to reduce a full 
infinite hierarchy to a single closed equation, which can be solved both numerically Q and analytically From the 
mathematical viewpoint, derivation of the BK hierarchy is just a method of reduction of the JIMWLK Hamiltonian 

, on some subspace of the functional space. In this paper, we try to generalize such a method into a wider class of initial 
J> ■ projectiles and consider a dipole that consists of a parton and an antiparton of arbitrary charge. This step allows 
| us to see the rich mathematical structure that arises in scattering amplitudes due to usage of the non-Abelian gauge 

■ group SU(N) in QCD. A physical application of our method is an explicit demonstration of intensive color duplication 
of Regge poles in a scattering amplitude. Higher Pomerons, which are associated with higher representations of the 
gauge group, have been observed recently [f|. They can arise in the analytical structure of the scattering amplitude 

' when the initial projectile is more complex than an ordinary dipole. In Ref. Q, where a gluonic dipole was studied, 
the higher Pomerons were considered. The method, which was used in is bounded with specific properties of 
adjoint representation. In this paper, we propose a formalism that allows us to equally study a dipole which consists 
of a parton and an antiparton of arbitrary charge. Mainly, we are interested in generalization of the BK hierarchy. 
The Pomerons can be analyzed in the weak-field limit of the BK equation. 

Usually, in applications the Balitsky-Fadin-Kuraev-Lipatov (BFKL) Pomeron Q and its corrections are widely 
exploited. It is believed that the BFKL equation explains thr initial fast growing of gluon density at small x. 
However, poles with a higher color charge potentially may influence the calculation of experimental quantities. From 
the theoretical viewpoint, a complete solution of evolution for an arbitrary projectile requires a formalism which can 
manipulate a contribution from any pole. 

An additional interesting question is about theoretical aspects of the origin of the BFKL equation, especially the 
question about its domain of applicability. Let us describe shortly the current known approaches to the BFKL 
equation: 

• In a classical approach to the BFKL equation a Pomeron is considered as a bound state of two Reggeized 
gluons. The disadvantage of this method is that it is not clear how to relate the BFKL equation with the 
full scattering amplitude of an arbitrary target and projectile. The bound state is constructed by hand, and 
there is no explicit algorithm of construction of the full amplitude. A so-called Reggeon field theory should be 
developed. 
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• We can use the dipolc model where the limit TV — > oo is assumed @. In this limit, the BFKL equation arises 
without problems as the limit of a small amplitude of dipole scattering. The disadvantage of this method is 
that the limit N — > oo is strongly unnatural and the rich algebraic structure associated with the realistic SU (3) 
gauge group is lost. 

• From the BK hierarchy we can derive the BFKL equation by applying the mean field approximation and the 
limit of weak dipole scattering amplitude. The mean field approximation reduces the infinite hierarchy to a 
single nonlinear BK equation. The disadvantages are that such an approximation has a limited domain of 
applicability and it is not clear how to generalize it to arbitrary initial projectiles. The advantage is only that 
target fields can be large without any restrictions. 

• The weak-field approximation is an assumption that target field distribution is concentrated near zero field. 
The method allows us to study complicated projectiles and more complicated poles such as the odderon We 
use this method in the current work in the process of linearization of a generalized BK equation. However, the 
general question here is about the domain of applicability of the weak-field approximation. Since the target 
field is weak, the target must be dilute. Hence, the dense-dilute picture is lost, and the JIMWLK equation is not 
applicable. In the dilute-dilute regime we must use another approach such as was proposed in However, 
in the current paper, we assume that the target field is small but is larger than the projectile field. So we use 
simultaneously both the JIMWLK equation and the weak field. 

This paper is organized as follows. In Sec. Owe develop the method of indexless transformation of color expressions, 
which is widely used in the paper. The method is simple and is based on decomposition of the Casimir operator on 
a tensor product of two representations. It allows us to significantly simplify subsequent calculations. As an example 
of the power of the developed method, in Sec. IHII we apply it to make a fast calculation of the nontrivial color factor 
which is a convolution of eight structure constants. In Sec. II VI we take a single parton with arbitrary color charge, and 
from the JIMWLK equation we derive the single parton evolution equation and its restriction to the fundamental case. 
We observe that this equation is a natural precursor of the common BK equation. In Sec. IVl from the single parton 
evolution equation we prove gluon Reggeization in an arbitrary color channel. In Sec. IVII we derive a generalization 
of the BK equation by taking a dipole which consists of a parton and an antiparton of arbitrary charge. Selecting 
fundamental representation, we easily reproduce the usual BK equation without any manipulation with color indexes. 
In Sec. IVIII by using weak-field approximation, we show how to obtain the BFKL equation. Section [Villi contains 
our conclusions. 



II. CASIMIR ON A TENSOR PRODUCT 



In many calculations in QCD we often need to transform a term like TJjMTJj, where R is some color representation, 
Tjj is a generator, and M is some matrix. Usually, such transformations are plagued by complicated index algebra. 
We develop the method of indexless decomposition that is based on the existence of the quadratic Casimir operator, 
which for irreducible representations is proportional to the unit operator. 

Consider a tensor product of two irreducible representations: A ® B. For the Lie group, generators in A ® B have 
the form 

T% 9B =id®T%+T%®id (1) 
The product A® B splits into a direct sum of irreducible representations: 

A ® B = V Q (2) 
Q 

where Q denotes representations and Vq corresponds to the invariant vector space. Generators also can be decomposed 
as 

t a®b = ^Z T Q ( 3 ) 
Q 

where Tq obey the usual reducibility conditions 

ImT% c V Q 



(4) 
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From definition JT]) the Casimir 1 operator for A (g> B is 

t a®b t a®b = C A + C B + ® (5) 
On the other side, from Q and (f4]) we have 

Q 

where Pq is a projection operator on subspace Vq. Index Q runs over all irreducible representations which belong to 
the decomposition of the tensor product. Projectors obey natural properties 

P 2 Q = p Q 

PqPr = 0; i£Q^R (7) 
Sp(P Q ) = D Q 

where Dq is the dimension of representation Q. Finally, we have 

2T% ®n = Y. °Q P Q - c a-Cb (8) 
Q 

Now consider representation R and any matrix M that acts in Vr. Under a gauge transformation matrix M transforms 
as M — > UMU' . Hence, it transforms as R®R : where R denotes complex conjugate representation (decomposition R® 
R = 1© other gives the natural correspondence between forms on Vr and vectors in Vg). Small gauge transformations 
in R have the form $ -> i>(l ~ iz a T a ). We can conclude that T| = -(Tr) t . Correspondingly, the term TgMTg can 
be viewed as 



TftMTft = (-Tg ® r|)M (9) 

The Casimir operator for complex conjugate representation is the same: = Cr. So with matrix notation we can 
write 

T a R MT a R = C R M -\Y, CqM q (10) 
1 Q 

where Mq = PqM. 

Equation (|10| is useful for transformation of numerous terms in QCD calculations. Its sufficient advantages are 
indcxless and simple generalization for arbitrary representations. Widely used in literature, 2 simple version of (|10p 
for fundamental representation is 

T ij T ks = T^StsSjk - J^SijSks (11) 

We shall use (fT0|) for derivation of a generalization of the BK equation for the dipole which consists of a parton and 
an antiparton of arbitrary charge. In the next section we shall demonstrate the method by fast calculation of a cube 
diagram for gluons. 



III. GLUON CUBIC DIAGRAM 



There are works where universal tools for multigluon color factor calculations were studied [ll|, [12J . See also [13j 
for computer friendly color flow decomposition. However, it would be suitable in partial situations to develop more 
special and simple method such as we have offered in the previous section. As an example, we calculate the nontrivial 
color factor in the gluonic Feynman corresponding diagram from Fig. [TJ which is a convolution of eight structure 
constants 



J _ yaia 2 a 3 ja 2 aia 5 ja 5 a 9 a e ja 3 a e a 7 jaja^a^ ja s a 9 a w ya4dloail jaianai2 



(12) 



1 In this paper we work only with the quadratic Casimir operator. 

2 It is called the "Fierz identity" 



FIG. 1: Cubic diagram. The vertex corresponds to structure constant f abc , and the edge is an index convolution. 



It was calculated in Ref. |14| by a birdtrack method which is unavoidably lengthy in this case. Of course, invariant 
(fT2"]) can be calculated using conventional properties of f abc and d abc symbols. Basic properties of these symbols can 
be found in the appendix to Refs. @, [H[. But this needs careful treatment of tensor expressions, especially the sign 
factor. Here we present the method of calculation in terms of SU(N) invariants. 
It is clear from Fig. [T]that the value of lfl~2|) can be expressed as 



I = Sp(Tg®Tg) A 

where we used the common fact that group structure constants are generators of the adjoint representation, 
following formula is a direct consequence of © 



T a R ®n = \Y J c Q p Q -c I . 



E 

Q 



C„ P 



(13) 



The 



(14) 



By inserting (fT4"|) into (fTB")) and using properties 0, we obtain 

I = Y.{\ C Q- C *) D Q 



Q 



(15) 



where Dq is the dimension of the representation Q. Now we need to know the properties of representations entered 
into the sum in (fT5|) . Decomposition of the tensor product of two adjoint representations of SU(N) group can be 
obtained using the Young diagrams. The Casimirs Cr can be calculated in various ways [l6|. Here we only collect 
and write in Table [J the required values taken from 
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1 N 2 - 1 


N 2 (N+3)(N-1) 
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2(N + 1) 
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2{N - 1) 
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10 

(JV 2 -4)(jy 2 



10 

-1) (JV 2 -4)(iV 2 



4 

2N 



4 
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TABLE I: Properties of representations entered into (|15[) . 



After some simple algebra we obtain the result 

N 2 

I = — {N 2 + 12)(N 2 - 1) (16) 
8 

There are two interesting consequences of ((T6|) . The first is the inapplicability of the large N c limit for the considered 
color factor in the SU(3) case due to factor iV 2 + 12. The second consequence is the S'[/(2)-scaling violation 3 . This 
scaling is the observation that in many types of diagrams the color factor for the SU (N) case is equal to the SU (2)-case 
factor multiplied by some integer power of (N c /2). 

The advantage of the presented calculations is the absence of any ugly tensor algebra. All relevant invariants for 
the SU (N) group can be collected once into compact tables and repeatedly used later. 



3 This clause is not about factor TV 2 — 1, which is just the dimension of the adjoint representation. 
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IV. SINGLE PARTON EVOLUTION EQUATION 

In many papers concerning the BK equation, authors usually start from the dipole scattering amplitude. However, 
we can also start from the single quark scattering amplitude. In this paper, we consider the case when a projectile 
consists of only one parton in arbitrary color representation R. The high energy evolution of parton wave functions is 
governed by the JIMWLK equation [l|, Q ■ Here we use the notation taken from Ref . [ijj • The JIMWLK equation is 

dS [ a ] rr, 6 lor I 



dY ' Sa 1 



zxy 



, v [-JZ(x)J$(y) - J a _{x)J a _{y) + 2V ha {z)j\(x)J a _{y)] 



_ {z-y)(z-x) 



ig Sa a (x, ±oo) 

+ CXJ 

V{x) = Pe — 

where Vb a (z) is the gluon scattering amplitude in an external field a a , T^ D are generators of the gauge group in 
the adjoint representation, and S[a] is the projectile scattering amplitude as a functional of target fields. In order 
to obtain an observable scattering amplitude, we need to perform an average over target fields with a corresponding 
weight functional. However, for theoretical purposes it is useful to use a nonaveraged functional S[a] with fixed target 
fields. Note that here we assume that the projectile is left-moving. 

Now we would like to find the action of the JIMWLK Hamiltonian on a single parton S matrix. The latter is given 

by 

S R (x) = Pe lg J a " {S ' x+)T * dx+ (18) 

where we assume that S(x) is a color matrix acting on the projectile color index. Functional derivatives can be easily 
evaluated: 

J$(x')S R (x) = 6(x' - x)T^S R (x) 

J1(x')S r (x) = 6(x' - x)Sr(x)T% (19) 

where are the SU (N) generators in representation R. Convolution with Vb a can be evaluated with the help of the 
following property which holds at any transverse point: 



After all evaluations we arrive at 

dS R (x) g 2 f 1 



dY (2tt) 3 J (z-x) 2 
Using decomposition PU|) . we have 

dS R (x) _ g 



S + T£S = V ab T b R (20) 



[2S R (z)T-S+(z)S R (x)T- - 2C R S R {x)} d 2 z (21) 



tt)3 J {z~x) 2 ^D R 



dY (2tt) 3 J (z-x) 2 

Q 



S R (z) (sUz)S R (x) 



d z z (22) 
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In the fundamental representation we can simplify Eq. (|22|) . We know that 3 ® 3 = 1 © 8. The Casimirs are C\ = 
and Cg = N. We have two projectors which obey Pi + P& = 1. A projector to an invariant state is Pi-M = Sp(M)/N. 
So for any matrix A we have 

A 8 = A-Sp(A)/N (23) 
By substituting this into (|22]) . we arrive at the first equation of " charged BK hierarchy" 

dS(x) _ Ng 2 f [S(z)S(x,z)-S(x)]d 2 z (24) 



dY (2tt) 3 J (z-x)' 
where 

S(x,z) = ±Sp(S+(z)S(x)) (25) 

In the last equation we can easily recognize a dipole scattering amplitude. Equation Q24[) has a very natural and clear 
physical meaning. A quark at transverse position x emits a gluon into the position z. The emitted gluon can be 
viewed as a quark-antiquark pair. The antiquark part of the gluon combines with the original quark into a dipole. 
Then the first term in (|24p can be viewed as multiple scattering of a dipole and quark component of an emitted gluon. 
The second term corresponds to the virtual correction due to the requirement of the overall probability conservation. 

It is not surprising that Eq. (|24[) is not closed. This means that functional S(x) does not form a complete space 
of solutions, and we must add corresponding equations for functional S(z)S(x, z) and so on. So there is an infinite 
hierarchy of equations. This hierarchy is just an attempt to reduce the JIMWLK Hamiltonian by some separable 
subspace of full functional space S[a\. If we take functional space V such as HV G V then the JIMWLK evolution 
can be reduced by V. The method of generation of hierarchy, like that started from ([Ml) , is the following. For given 
starting space V) such as the space of S(x) in (|24[). we generate space V = ® n H n Vo. It is clear that HV G V, so the 
evolution can be reduced by V. 

There are two approximate methods for closing an infinite hierarchy to a finite number of equations. The first method 
is the mean field approximation where the target average leads to something like (S(z)S(x, z)) — > (S(z))(S(x, z)). 
The functional dependence on a a (x) is removed, and we deal only with ordinary functions. The second method is the 
weak scattering approximation where we reduce full functional space by considering functionals only on small fields 
a a {x). Since the physical functional S[a] must obey S[0] = 1, at small a a (x) the functional 1 — S[a] is small, too. 
We apply this method in the next section. 



V. GLUON REGGEIZATION 



Now we want to switch to the weak-field limit where it is assumed that the target fields are small. In this section, 
representation R is assumed for symbols S and M. In the weak scattering limit we define scattering amplitude M(x) 
as 

S(x) = 1 - M(x) (26) 

In this limit M(x) < 1, so we can keep in (1211) only linear over M terms and can use identity M ] = —M. By 
expanding (|2~Tj) , we obtain 

dM{x) _ g 2 [2C R (M(z)-M(x))-T^(M(z)-M(x))T^]d' z z (27) 



dY (2tt) 3 J (z-x) 
By substituting the expression @ into (|2T|) and making use of ([5]), we finally arrive at 

dM Q (x) _ C Q g 2 f 1 



dY (2tt) 3 J (z-x) 



M Q (z)~M Q (x)]d 2 z (28) 



for each irreducible representation Q which contributes to a tensor product R®R. The matrixes Mq are Mq = PqM 
as it was defined in (|10[) . 

Equation (|28|) provides the closed expression for the evolution of the scattering amplitude of the projectile in the 
given representation R in the channel with the given color exchange Q. To find a complete solution, we need to 
decompose the initial conditions of M to a sum of irreducible representations and solve for each component Eq. 
Note that the singlet representation has zero Casimir, and it gives a constant solution of 
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Equation (I28|) can be easily solved in the momentum space. By converting f|28[) to the momentum space, we obtain 

dM Q (k) 



dY 



LJ Q (k)M Q [k) 



The solution can be easily obtained as 



M Q (k, Y) = M ( °\k)e^^ Y (30) 



where M { ° ] = PqM^ are initial conditions of evolution. This solution has Regge form M ~ s a ^ as it should. If Q 
equals the adjoint representation, then the known expression of a gluon pole trajectory has the form [l8| 

= -ITT I — ^ d "P ( 31 ) 

Integrals in (|2"i?|) and (f3"Tj) can be evaluated and they give equivalent answers 

where /i is an infrared regulator. 

There is one important subtlety here. Though in the weak-field limit we require M <C 1, individual components 
of M may have distinct order in comparison with each other. Indeed, let us recall that M = 1 — exp(ia a T a ). Since 
there are two obvious identities (T a )§ = T a and (T a )^§ = (index 8 denotes here the adjoint representation that is 
constructed from generators), we have 

M# = 0(M|) 
M ~ M 8 

This means that the main contribution in (|2"5)) comes from the adjoint representation with Cs = N. The other Mq has 
order at least (M§) 2 , which is negligible in the weak field limit. However, we can study Eq. 1|28[) beyond the weak-field 
limit, too. It may have sense as a term of formal power expansion of the full evolution equation (J22J)- This expansion 
may be useful for construction of something like the Reggeon diagram technique. Also, higher representations Q may 
be relevant when the target has specific field correlators which allow one to set (MgAfg) ~ and similar for higher 
powers. 

The evolution equation (|29p corresponds to a single moving pole in a complex angular momentum plane. The pole 
trajectory is given by solution (|32[) . So we can see the so-called reggcization phenomenon. However, we have found 
more than one pole - one pole for each irreducible representation Q. This is a really remarkable result since there is 
a transfer of gauge group algebraic structures to the analytic structure of the scattering amplitude. Conversely, from 
the analytic structure of a scatteringamplitude we can enumerate irreducible representations. A similar situation is 
observed in the Pomeron trajectory [6||. Since the Pomeron can be viewed as a bound state of two Reggeized gluons, 
we naturally conclude that there is an intensive color duplication of poles. Poles with higher Q can be relevant 
where there are many native partons (quarks or gluons) in a projectile. Several native partons can form a higher 
representation via the tensor product of their color spaces. The significant feature of the BFKL Pomeron case is that 
there are many poles even without color duplication. This happens due to strong degeneration of the eigenvalues of 
the spectrum of the BFKL operator. 

VI. GENERALIZED BK EQUATION 

Now we consider the case where the projectile is a color dipole built from two partons in representations R and R. 
We want to study evolutions of the following scattering functional: 

S=^-Sp(uHy)U(x)) (34) 
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where U is the Wilson line of representation R. Since the JIMWLK operator H in (fl~7| is a differential operator of 
second order, in calculation of HS we can apply the Leibnitz rule, which gives four terms which are equal to each 
other but with different spatial kernels. We show the term where all J a act on the U(x) in (f34|) . It is 

±-Sp{UlHU x ) = J 0^ Sp {-UlT a R T a R U x - VlV x T%TH + 2V ba UlT b R U x T a R ) (35) 

The calculation of the other three terms is very similar. By using T R T R = C R and the identity (|20[). which is read as 

V ba (z)T b = U{z)T a U\z) (36) 

and adding three remaining terms, we obtain 

% = HS= J^fj M '*»^ Sp WU z T a R UtU x T^ - 2C R UlU x ] (37) 

z 

where M zxy is the well-known dipole kernel 



■M ZX y — K zxx -\- K Z yy zxy -^zyx ~~ ( \2( \2 (^^) 



(x - yf 

(z — x) 2 (z — y) 2 

The key point of our method is the usage of decomposition (|10p . which in the current context has the form 

T a R UtU x T a R = C R U\U X \ C Q ( U 1 U -)q (39) 

Q 

As usual, Q runs over irreducible representations which contribute to the tensor product R® R. By substituting (|39j) 
into p7p . we finally arrive at 

§ = j£pf M -E [-UlU 2 (UlU x ) Q ] (40) 

When R equals the fundamental representation, we can easily obtain the usual BK equation. Using Eq. (f2"3")l . we 
obtain a simplified version of (|40| : 

f M zxy [S(y,z)S(z,x)-S(y,x)} (41) 



dY ~ (2tt) 3 J zxv 

z 

Equation (|40[) is the first equation of a complicated nonlinear hierarchy which has an additional complexity level in 
comparison with the usual BK hierarchy. This happens due to the presence of the sum over various Q projections. 
When we further calculate HSy Z S zx we obtain new functionals which are various combinations of the Wilson lines, 
and when we perform a corresponding operation on right-hand side of (|40|) we obtain such functionals with two various 
Q projections. It should be stressed that linear and nonlinear terms in (|4"Tj) are single whole which is just a projector 
on the adjoint component of the tensor product. This presents a contrast to the common description where the linear 
term is explicated by so-called virtual corrections. 



VII. BFKL EQUATION 

At first, let us show shortly how to obtain the BFKL equation from (|41[) by using the mean field approximation 
as mentioned in the introduction. We should not forget about target averaging in (|4ip . The approximation gives 
(S(y, z)S(z, x)) = (S(y, z))(S(z,x)). Next, we define (S) = 1 — N. By taking the limit N <C 1, we obtain the usual 
BFKL equation. It should be stressed that claim JV <C 1 does not assume a small target field. Unfortunately, it is 
not clear how to generalize the considered method to a case of more complicated projectiles. 

In order to obtain the generalized BFKL equation, we consider the weak scattering limit of (|40[) . Let U = 1 — M. 
Since U E SU(N), we have 



JU 1 + M = MM** 



(42) 
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Since Mi = Sp(M)/ Dr and C\ ~ 0, wc have the following useful relation for any M and Q: 

C Q Sp(M Q )/D R = Cq(Mq)i = (43) 
and there is a similar fact about the unit matrix 

CqIq = (44) 
By expanding up to second order the matrixes in the right-hand side of (|40|) over M , we have 

XJp z (Up x ) Q = 1 Q - Mtl Q - M z 1q - (Mt)Q - (M x )q+ 

+M}M z l Q + Mt(M]) Q + Mt(M x ) Q + M z {Mt) Q + M Z {M X ) Q + (M]M X ) Q 



Next, using (J42J) , (|43J) , and (|44J) we obtain 
diV(y,a;) <? 2 



M zxy V C Q [JV q (i/, z) + AT Q (z, x) - Nq(u, x) - N Q {z, z)} (46) 



where we have defined 

N Q (y,x) = 1 - 5 (y >a r) = 1 - ^-Sp(U^y)U Q (x)) = ^-Sp(M^y)l Q + M Q (x) - M\y)M Q (x)) (47) 

It can be easily checked that the linear over M terms cancel each other in (|4"6l) . 

It is instructive to obtain the original BFKL equation from (|4^|) . Note that before this point we used only the 
assumption M <C 1. However, as was shown in (|33l) of Sec. [V] in the weak-field limit only the M§ component is 
relevant. The other Mq has order at least 0(M|). Hence, the sum over Q in (|46|) is reduced to one term with Q = 8. 
In additional, properties (|33|) allow one to set Ng(y,x) = N(y,x). So we obtain the original BFKL equation 



dN(y,x) _ g N t 



dY (2tt) 3 



M zxy [N(y,z) + N(z,x)-N(y,x)} (48) 



where we used N(z, z) = 0, since U Z U\ = 1 exactly. 

Beyond the weak-field limit, from a power expansion viewpoint, Eq. (|46[) inevitably contains the term Nq{z,z). 
Note that a similar term is already known. In Ref. Q it arises in the evolution equation for the 2-point Green's 
function in the weak-field regime. The functional Nq{z, z) depends on the field a a {z) and can be viewed as a function 
of one variable on a group manifold. The mapping a a — > SU(N) is cxp(io! a T a ). This function on a group manifold is 
intrinsic, like characters, in the sense that it can be constructed directly from the definition of a group. The system of 
equations (|4C5|) is not closed because there are many unknown variables Nq in the right-hand side. This means that 
there are exchanges of many different poles in the dipole scattering amplitude. To find the color diagonal version of 
([4"6")h we must consider the action of the JIMWLK operator H on a two parton matrix element <g> U x with arbitrary 
color indexes and perform a diagonalization of color structure. In particular, we can show that the system of equations 
which is obtained from evaluation of HSp(U^(U x )q) in the weak scattering approximation is linear and closed. This 
means that the evolution equation for Nq(jj,x) contains only other NQi(y,x). Unfortunately, this system is intricate 
in comparison with the original BFKL equation. Sec Ref. [6j for the eigenvalue problem in the adjoint case and Ref. 
[HI for a detailed study of multigluon states. 



VIII. CONCLUSION 



In this paper, we have studied one relation of gauge group algebraic structures to the high energy QCD evolution. 
The classification of classical Lie groups and its irreducible representations is the classical mathematical result which 
is widely used in modern theoretical physics. We are convinced that a full solution of high energy QCD evolution 
must be naturally related to the rich mathematical structures on the gauge group. In fact, many known investigations 
ignore such structures. We have seen that the form (|4"Tj) of the usual BK equation is not natural and the more 
natural form is Eq. ((40]) . Moreover, it is clear that the initial choice of a color dipole is not natural either. More 
simple and transparent equations emerge when we start evolution from a single parton. Geometrically, such choices of 
initial states are equivalent to selecting finite polynoms of matrix elements of irreducible representations as an initial 
condition of functional S [a]. It is clear that during evolution the power of such polynoms is growing. So it is difficult 
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to find appropriate analytic solutions for all Y . Currently, our road map is to work on a functional level [l9| . Powerful 
mathematical techniques can be used if we view S[a] as functions on the group manifold. 

One can think that higher representations are irrelevant to the real world, where we deal only with quarks and 
gluons. However, if a projectile has many partons and we use coarse transverse resolution, then few partons via tensor 
product can form a new effective parton, which belongs to the representation of a higher color charge. This fact was 
widely used in Ref. (l6| . where it was shown that in the dense case the charges with higher Casimir become dominant. 
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